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Screened Vortex Lattice Model with Disorder 

C. Giardina 1 - 2 , N. V. Priezjev 2 , J. M. Kosterlitz 2 
Dipartimento di Matematica, Universitd di Bologna, 40127 Bologna, Italy and 
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The three dimensional XY model with quenched random disorder and finite screening is studied. 
We argue that the system scales to model with A ~ ~ T and the resulting effective model is 
studied numerically by defect energy scaling. In zero external field we find that there exists a true 
superconducting phase with a stiffness exponent 8 ~ +1.0 for weak disorder. For low magnetic 
field and weak disorder, there is also a superconducting phase with 9 ~ +1.0 which we conjecture 
is a Bragg glass. For larger disorder or applied field, there is a non superconducting phase with 
6 ~ —1.0. We estimate the critical external field whose value is consistent with experiment. 
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The phase diagram of type II superconductors in an 
external field has been the subject of intense theoretical 
and experimental investigation |Q. After the discovery 
of high-T c materials, the role of both thermal and disor- 
der induced fluctuations has been reconsidered, revealing 
many new interesting phenomena. In clean systems, it 
was realized [0 that, with increasing temperature, the 
Abrikosov lattice melts into a vortex liquid via a ther- 
mally induced first-order transition. Experiments per- 
formed on thermodynamic quantities such as magneti- 
zation || and specific heat B confirmed the first order 
nature of the melting transition in YBCO and BSSCO 
materials. More recent studies of the clean unscreened 
system in an external field || have shown that the low 
temperature vortex lattice melts at T = Tm to a liq- 
uid of rigid flux lines and at Tr, > Tm the lines become 
entangled and vortex loops proliferate. Very recently, 
there have been studies || 0, || of a random 3D XY 
model wthout screening. These claim to study the glass 
transition in a superconductor but do not have the vital 
screening and obtain results which are inconsistent with 
each other so their relevance to real systems is unclear. 

In this Letter, we study the stability of superconduc- 
tivity in the presence of point disorder and an applied 
magnetic field by computing numerically the stiffness ex- 
ponent 9. We argue that the effective screening length 
A and temperature T scale to zero at very long length 
scales so that a model with A ~ ~ T, which is amenable 
to simulation, is a physically reasonable model of a dis- 
ordered superconductor at low T. Using this effective 
model, we can identify the field and disorder driven tran- 
sition H from a superconducting to non-superconducting 
phase as the field or disorder is increased and conjecture 
that the former is a Bragg glass |L0]]. The nature of the 
large disorder phase is still controversial: a viscous non- 
superconducting pinned vortex liquid [ pd| or a supercon- 
ducting vortex glass |TJ] . 

The ingredients necessary to describe a typical high-T c 
superconductor in a field are (i) pinning of flux or vor- 
tex lines and loops by random point impurities and (ii) 
weak screening of the interactions between vortices. We 



argue that a model containing these essential ingredients 
is a three dimensional XY model on a simple cubic lat- 
tice with quenched random phase shifts. In the vortex 
representation, the Hamiltonian is [l5[ 



(1) 



We may ignore boundary terms which, at least in the best 
twist approach, vanish by a proper choice of global twists 
p6| . The dynamical variables are the integer valued vor- 
ticities J \ on the links of the dual lattice and subject to 
the local constraint (V ■ J)i = at every site i. The 
bi arc quenched random fluxes on the dual lattice which 
are obtained from the circulation of the quenched vector 
potential A and by adding a uniform external field h in 
the z direction 



b t = -L[V x A]i + hz 



(2) 



Here h = Bciq/Qq is the mean flux per elementary pla- 
quette normal to the applied field Bi, a® is the underly- 
ing lattice spacing and <I>o = 2 x 10~ 7 Gauss-cm 2 is the 
flux quantum. The vector potential A^i with ji = x,y, z 
is independently uniformly distributed A^i € [0, 1~kol) 
with < a < 1 and is defined on the bonds of the origi- 
nal lattice. The disorder strength a interpolates between 
two well known limits, the clean system, a — 0, and 
the maximally disordered, a = 1, gauge glass. By con- 
struction, the fields bi satisfy the divergenceless condition 
(V • b)i = on every site. G(i,j) is the screened lattice 
Green's function with dimensionless screening length A 
in units of eto 



G(i,j) = 



£3 



E 



exp[ik ■ (rj - rj)] 
2£Jl-cosfc A1 ) + A- 



(3) 



= (xi,yi, Zi) is the i-th site on the dual lattice 
n — lnn^jL, with /i = x,y, z and rin = (!,••■, L). 



where r 
and fc 

The vortex system defined by Eqs. (|l|), (g) and (||) is a 
vortex lattice of spacing ft, -1 / 2 for a — in the absence 



2 



of disorder. Turning on disorder distorts the lattice and 
pins it at a favorable position in the underlying substrate 
lattice, as is clear from Eq. (H), but does not destroy 
the underlying periodicity of the vortex structure [ p)[ . 
Real cuprate superconductors are very anisotropic with 
r = K.IK,y ~ 5 for YBCO and T « 100 for BSSCO 
jl| where we have chosen the z direction as the c axis 
and x, y directions as the a,b axes. We note that when 
X u <C 1, the vortex- vortex interactions become isotropic 
||17| and, in this limit, our isotropic model is realistic. 

The final and most important part of our argument is 
to justify the validity of the strong screening limit for an 
extreme type II superconductor where the bare screening 
length oo > A 3> 1 is large but finite. ^From Eq. (||), it is 
clear that A scales like a length so that scaling ao — + e'ao 



induces the scaling A 



,-z 



A. Thus, at long length scales 



which are relevant to the weakly disordered system of 
interest, the effective screening length A <C 1. In this 



strong screening limit, G(i,j) 
yields 



(27rA) 2 5y+0(A 4 ) which 



H 



(2ttA) : 



Jew 



b I ) 2 + 0(A 2 



(4) 



This local form of Eq. ([!]) on a cubic lattice can be re- 
garded as a model for the long length scale properties 
of a disordered superconductor and can be studied nu- 
merically by very efficient combinatorial optimization al- 
gorithms jl8|, Q on large systems. Any nonlocal terms 
such as vortex - vortex interactions or boundary terms 
|l6| , |20|] in Eq. (|l|) render such algorithms useless. Of 
course, screening is weak, A > 1, in a real type II su- 
perconductor and one may question the relevance of a 
system described by Eq. (Q). Since our interest is in the 
superconductivity and in the vortex lattice structure as 
a secondary effect, we argue that Eq. (Q) is an adequate 
description of a superconductor in an external field with 
any finite screening. Numerical simulations jl6[ Q 
of Eq. (0) indicate that the exponent 9 has the same 
value for any A < oo. We find this slightly surprising as 
decreasing A transforms a type II into a type I supercon- 
ductor and the independence of 9 on A seems to imply 
that these do not differ in any essential way at small 
T. Additional justification for the relevance of the model 
of Eq. (^) at T = is that T is an irrelevant variable 
pOf scaling to zero in the superconducting region. This 
gives some justification for the computationally accessi- 
ble model of Eq. (Q) and for our scaling arguments. 

To investigate whether a transition occurs, we use de- 
fect energy scaling. In this approach, one computes the 
energy AE(L) of a defect in a system of linear size L and 
fits to the scaling ansatz 



(AE(L))~L 6 



(5) 



where 9 is the stiffness exponent and (■ • •} denotes an 
average over disorder. The sign of 9 distinguishes two 



regimes: if 9 > 0, inserting a defect costs an infinite en- 
ergy in the thermodynamic limit and the system will be 
ordered at sufficiently small finite T . Conversely, if 9 < 0, 
large domains cost little energy and, at any T > 0, su- 
perconductivity will be destroyed. To calculate the de- 
fect energy we employ the method proposed by Kisker 
and Rieger |2^] , who restated the problem of finding the 
ground state for Hamiltonian (ji|) in terms of a minimum- 
cost- flow problem [fl8| , where the cost functions are pre- 
cisely given by Ci(Ji) = (Ji — bi) 2 /2. This method makes 
use of the successive shortest path algorithm (SSPA) (n]] 
to find the ground state configuration {J } for each re- 
alization of disorder. The global flux / associated with 
this configuration is given by 



(6) 



The elementary low energy excitation configuration 
{J 1 } is obtained by gradually decreasing all costs in, 
say, the z direction until the global flux f z jumps by 
one, f z — y f z + 1. The lowest energy excitation will 
be a global vortex loop encircling the 3D torus in the 
z direction. The defect energy is then obtained by 
AE = E^J 1 }) — E({J }). A more conventional way 
of determining the defect energy is to calculate the en- 
ergy difference between periodic and antiperiodic bound- 
ary conditions, which amounts to adding a global twist 
of 7r along one spatial direction jll| . In our case we can 
ignore the boundary terms because the global loop corre- 
sponds to a twist of 27r, which has no effect as the original 
Hamiltonian is invariant under a discrete gauge transfor- 
mation modulo 2ir. Remarkably large system sizes can be 
treated by applying the SSPA, while conventional meth- 
ods such as repeated quenching or simulated annealing 
are much less efficient. In this work we study L x L x L 
systems with L < 40 for different values of a in the range 
< a < 1 and magnetic field in the range Q < h < 0.25. 
The number of realizations of the random bonds varies 
from 500 for 40 3 systems up to 10 4 for the smallest ones. 
There is no upper limit on L. 

The zero field case has been studied in detail |^3| and 
we repeated the simulation to check our algorithm. The 
results are summarized in Fig. ([!]) for several values of 
a and they are essentially identical to those of Ref . [^3| . 
For weak disorder, a < a c ~ 0.495, we obtain 9 ~ +1.0 
indicating a superconducting phase for T > 0. For strong 
disorder, a > a c , the stiffness exponent 9 ~ —1.0 indi- 
cating a non superconducting phase for T > 0. 

In an applied field type II superconductors have a fixed 
density of vortex lines which, with weak disorder, form a 
distorted Abrikosov lattice or Bragg glass |l(| at small T. 
Increasing the field effectively increases the disorder and 
the Bragg glass phase transforms into an entangled ar- 
ray of vortex lines with no remnant of translational order. 
The model of Eq. p) contains the standard description 
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FIG. 1: Size L dependence of domain wall energy {AE) in zero 
external field (log-log plot). The legend shows the magnitude 
of the disorder strength. Solid lines are guides for the eyes, 
dashed lines with slopes ±1 are drawn for reference. 



of a vortex lattice as an elastic periodic solid in a random 
pinning potential when A > 0. The essential supercon- 
ducting properties seem to be unchanged as A — > but 
the interaction of the flux lines determining the periodic 
lattice structure vanishes when A = 0. Our strategy is 
to impose the required periodic vortex lattice structure 
by imposing a periodic array of columns of bonds along 
z which favor flux lines on these columns. We make the 
plausible conjecture that this procedure retains the es- 
sential physics of a disordered superconductor at T ~ 
and it is the only system which can be simulated with 
large size L. To compute the defect energy AE{L) it is 
essential to impose periodic boundary conditions in all 
directions and imposing a fixed number of flux lines by 
the source/target method [[l8| is not applicable. 

We implement the periodic lattice potential by first 
adding a uniform h — Ba^/^o to all i bonds. To en- 
sure that the vortex density does correspond to this h, 
we introduce a potential A on the bonds of the expected 
vortex lattice by h — ► h + A in Eq. (0) . The costs of the 
flux lines are decreased by increasing A until the match- 
ing is satisfied with one flux line per favorable bond. We 
thus obtain a ground state of energy Eq{L) with the re- 
quired number N = hi? of flux lines. The defect energy 
Ed(L) is obtained in the standard way (2^, Q by reduc- 
ing costs until an extra flux loop is induced. This is the 
simplest way of introducing the essential periodic part of 
the random pinning potential on the flux lines jTo| , It 
is interesting to note that our model is essentially that 
of many elastic lines in a periodic potential J2j| so the 
transition at a c (h) is a transition from a superconduct- 
ing state to a state which is both rough and entangled. 
This is in contrast to a superconductor in a field with no 
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FIG. 2: L dependence of domain wall energy (AE) for finite 
external field h = 0.25. Varying the disorder strength (see 
legend) changes the sign of the stiffness exponent 8. Dashed 
lines have slopes ±1, while solid lines are guides for the eye. 



disorder where the vortex lattice first loses translational 
order and becomes a liquid of rigid lines and these lines 
become entangled at larger field or temperature [EJ. 

In Fig. (Q) we show the behavior of the defect en- 
ergy with system size L for a fixed value of external field 
h = 1/4 for different disorder strengths a. It is more 
convenient to fix the field and vary abecause the pe- 
riod of the vortex lattice must be commensurate with 
the system size L which allows only restricted values 
of h = NL~ 2 with N, L integers. At small disorder, 
a < 0.40, we observe a positive stiffness exponent, which 
asymptotically tends to 6 = 1. Increasing the disorder, 
a > 0.40, we find that the defect energy decreases with L 
and for large L, (AE(L)) ~ L^ 1 . There is a critical value 
a c (h — 0.25) ~ 0.40 separating an ordered from a disor- 
dered phase. Fig. (||) shows typical ground state config- 
urations of the system. Below the critical disorder, Fig. 
(|) left, the lowest energy configuration forms an almost 
perfect vortex lattice, while above the critical disorder, 
Fig. (H]) right, the lines are rough and entangled ]25| , as 
expected for a phase with proliferation of dislocations. A 
similar analysis for different values of h shows that a c (h) 
decreases monotonically as h increases as shown in Fig. 
(^) . We estimate the critical value of the external field as 
h c = 0(1) where h = Bo^/^q is the mean flux per pla- 
quette. Here B is the actual field, ao is the lattice spacing 
and <f>o = 2 x 10 -7 Gauss-cm 2 is the flux quantum. In 
the strong screening limit, ao ~ A, whose typical value 
in high-T c superconductors is A ~ 10~ 5 cm. This yields 
an estimate for the critical field B c ~ 10 3 Gauss. This 
is to be compared with the typical experimental value of 
B c ~ 500 Gauss in BSSCO §. However, this should also 
hold for YBCO which has a comparable A at T = ffl 
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FIG. 3: Ground states for h = 0.25 in the ordered phase, 
a = 0.10 (left) and in the disordered phase, a = 0.48 (right). 



0.55 



0.50 



a " 0.45 



0.40 



0.35 



0.00 




FIG. 4: Critical disorder strength a c (h) for < h < 0.25. 



and our estimate of B c depends on A only. This seems to 
disagree with experimental estimates of B c ~ 10 4 Gauss 
fjfjf but the measurements are for T > 50°K. 

In this Letter we have studied the strongly screened 
vortex glass in the presence of disorder and, for the first 
time, have successfully implemented a defect energy scal- 
ing study of the stability of superconductivity as indi- 
cated by the stiffness exponent 9 in an external field us- 
ing periodic boundary conditions. It would be interest- 
ing to apply this method to probe directly the behavior 
of the model allowing for dislocations in the vortex lat- 
tice by applying appropriate boundary conditions p7| . It 
would also be interesting to treat a more realistic model 
with finite screening and to show more convincingly that 
the T = A = limit studied in this Letter is the stable 
renormalization group fixed point of a disordered super- 
conductor at low temperature. 
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